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Abstract 

A low energy effective field theory model for QCD witti a scalar color octet field is 
discussed. The model relates the gluon mass, the constituent quark masses and the 
quark condensate. The gluon mass comes about \fNc^QCD with the quark conden- 
sate being proportional to the gluon mass squared. The model suggests that the 
restoration of chiral symmetry and the deconfinement transition occur at the same 
temperature and that, near the transition, the critical exponent for the condensate 
is twice the gluon mass one. The model also favors the decoupling like solution for 
the gluon propagator. 
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The use of effective field theories is a common practice in Physics (see e.g. 
[1]). There are many examples in the literature of such type of approach for 
low energy hadronic physics from chiral perturbation theory (see e.g. [2]), to 
quark models with contact interactions and quark models including pions and 
meson fields (see e.g. [3]). In this work, we discuss an effective model, invari- 
ant under local SU{'i) color, for the non-perturbative regime of QCD which 
uses as degrees of freedom the gluon, the quark and a color octet scalar field 
(j)'^. The scalar field can be viewed as multi-gluon and multi-quark excitation 
of the vacuum, i.e. it should be interpreted as a collective excitation of the 
QCD vacuum. Given that the scalar field is not a color singlet, 0" does not con- 
tribute directly to the S-matrix. The scalar field resumes the non-perturbative 
properties of QCD and there is no need to specify in detail its dynamics. The 
0° contribution to the quark and gluon dynamics comes via the color singlet 
condensate which, according to our estimates, turns out to be of the 

order of Aqcd- The scalar field generate a gluon mass rrig, shifts the current 
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quark mass and contributes to the quark condensate (gg). It turns out that the 
shift in the quark mass and the quark condensate are all proportional to the 
effective gluon mass. In this sense, the model links the presence of an effective 
gluon mass with chiral symmetry breaking (xSB). 

The studies of the phase diagram of QCD [1] suggest that the restoration 
of chiral symmetry and the deconfinement phase transition occur at similar 
values for the temperature and chemical potential, at least for small enough 
values of the chemical potential /x < 300 MeV. In the model under discussion 
either the gluon gets a dynamically generated mass, quarks are confined and 
chiral symmetry is broken or the gluon mass vanish, quarks become deconfined 
and chiral symmetry is restored. The model predicts that the restoration of 
chiral symmetry and deconfinement transition occur simultaneously at low 
baryon density, in close agreement with the study of the QCD phase-diagram. 

Some of the results obtained here can be tested comparing the predictions with 
non-perturbative solutions of QCD, namely with the results from Schwinger- 
Dyson equations (SDE) and lattice simulations. For pure gauge, the recent 
solutions of the SDE can be classified into two classes. The scaling solution 
[S] which predicts a vanishing gluon propagator at zero momenta and has no 
connection with a dynamical generated gluon mass. The decoupling solution 
PITfS] which describes the gluon as a massive particle in the low energy regime, 
with the non- vanishing mass being related with gluon confinement. Other 
authors have also contributed to the understanding of the solutions of the SDE, 
see [H] . The two solutions agree in the ultraviolet region but differ significantly 
in the infrared region, with the decoupling solution being closer to the results 
of lattice QCD simulations. 

Lattice QCD simulations [lU] also suggest that the gluon behaves as a mas- 
sive particle in the low energy regime and becomes massless in the high energy 
regime. A non- vanishing gluon mass is welcome to regularize infrared diver- 
gences and to solve unitarity problems. Furthermore, a dynamically generated 
gluon mass is related with the presence of a (A^) gluon condensate. A non- 
vanishing gluon condensate is a clear sign of the non-perturbative sector of 
QCD - see, for example, [11] and references therein. Besides the theoretical 
support for a gluon mass, diffractive phenomena ^2] and inclusive radiative 
decays of J/ip and T [13] also suggest a massive gluon. The precise value for 
the gluon mass depends on the definition used to compute it. Lattice QCD 
simulations and Schwinger-Dyson equations suggest an infrared gluon mass 
of rrig ~ 600 MeV [6|10fl4j . Phenomenology points towards a gluon mass be- 
tween 0.5 GeV and 1.2 GeV, depending on how the mass is defined; see table 
15 in [15]. 

The effective model predicts a constituent light quark mass M oc m^. This 
result can be tested looking at the ratio M/m"^, as a function of momenta, using 
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the decoupling solution of the SDE. The comparison was performed using the 
results of [8], which rely on quenched lattice gluon and ghost propagators 
to solve the quark gap equation in the chiral limit. The decoupling solution 
depends on the ansatz used for the quark-gluon vertex. In [£| the authors 
explore two different definitions for the quark-gluon vertex. It turns out that, 
in both cases, the ratio of the constituent quark mass over the effective gluon 
mass squared is compatible with an essentially constant M/ rri^ in the infrared 
region. In this sense, the predictions of our model get some support from the 
non-perturbative solution of the QCD Schwinger-Dyson equations. 

In the formulation of the effective model we assume flavor independence to 
reduce the number of independent parameters. Flavor independence of color 
interaction has phenomenological implications that have been checked exper- 
imentally at SLAC \TE\, measuring the strong coupling constant as{Mz) in 
Z° — bb{g), — )■ cc{g) and QiQiid), where / = u,d,s. It was found 

that, within errors, as{Mz) is flavor independent. Furthermore, flavor inde- 
pendence of strong interactions was also explored on the theoretical side to 
understand hadronic phenomena [T7] . 

The effective model adds to the usual QCD Lagrangian new operators which 
are, from the point of view of the quark operators, either color singlet or octet. 
There is a unique quark singlet operator, which mimics the so-called ^Pq model 
[18] used to describe OZI-allowed mesonic strong decays. Recall that the ^Pq 
model assumes that a pair of quark-ant iquark is created out from the vacuum, 
with vacuum quantum numbers, to create new mesons by the recombination 
of the quarks in the initial state and the "vacuum quark-antiquark" pair. 

Effective Degrees of Freedom. Quantum Chromodynamics can accommodate 
multiquark and pure gluon states. From the point of view of an effective theory, 
these multiparticle states can be viewed as composite fields. The success of 
lattice QCD in the quenched approximation can be seen as an indication that 
the non-perturbative QCD dynamics is mainly in the gluon sector. At the 
level of the Lagrangian, it is precisely the gluon sector which contains terms 
of third and fourth order in the gluon field A^. We will assume that the non- 
perturbative physics is resumed in the simplest tensor structures, i.e. scalar or 
pseudoscalar fields, whose main contribution comes from gluonic operators. 

For global color transformations belongs to the octet representation. Given 
that 8 O 8 = 1 © 8 © 8 © 10 © To © 27, the simplest two gluon color operator 
is a color singlet. Lattice QCD simulations in the quenched approximation 
[19] predicts the mass of the low-lying glueballs to be 1710 MeV for quantum 
numbers J^^ = 0++, 2560 MeV for J^^ = 0"+ and 4780 for J^^ = 0+". We 
aim to build an effective theory describing hadronic phenomena below the 1 
GeV energy scale. From the values for the glueball masses, we expect that 
the scalar and pseudoscalar glueballs to play a secondary role and we will not 
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include color singlet scalar fields in the effective theory. 



From the gluon field itself one cannot build a color octet under local transfor- 
mations. Instead, one can use the non-abelian Maxwell tensor which, for 
local color transformations, belongs to the adjoint representation of SU{3). 
The two possible color octet fields with different Lorentz indices are 

•^fiwqc = fabcF^j^F^^ and S^^^iC, = '^"■^c F^^^F^^ , (1) 

where fate (dabc) is the antisymmetric (symmetry) structure constants of SU (3); 
see, for example, for definitions and properties. The color transformation 
properties are derived from fabc and date Jacobi identities. Note that in ([1]) 
one can replace the Maxwell tensor by its dual without changing its trans- 
formations properties. The composite fields A and S were defined to have 
the mass dimensions usually associated with a bosonic field. This requires the 
introduction of the mass scale A. 

The composite operators A and S generate multi-gluon states from the QCD 
vacuum. The contraction of Lorentz indices for the product of two Maxwell 
tensors gives rise either to a scalar field F'^^^F^^'^ or a pseudoscalar field 
F^T^F^^^ ■ Of course, besides the scalar and pseudoscalar composite fields, other 
Lorentz tensors can be built from ([T]). However, we are assuming that the 
scalars or pseudoscalars capture the essential of the non-perturbative glue 
physics and more complex configurations represented with higher rank tensor 
operators will be disregarded. The spin zero composite fields belonging to the 
color adjoint representation built from equation ([1]) are 

^dabcF'^,F''f'\ if r is a scalar, 
■p- dabc F^j^F'^^^ , if 0" is a pseudoscalar field. 

An operator with the same quantum numbers can be built using only quark 
operators. Therefore, one can write a scalar color octet composite field as 

r = ^46cF;.f^'^'^ + ^[gt'^g] + (3) 

where a and a' are constants. The dots in ([3]) stand for other possible operators 
including multi-gluon, quarks and ghost fields. If 0" is a pseudoscalar field, in 
the first term of ([3]) one should replace one of the non-abelian Maxwell tensors 
by its dual and in the fermion a 75 should be included. 

The definition ([3]) allows a naive estimation of the relative contribution of 
gluon and quarks to cp"' via the ratio TZ = {F^) / A{q q) , where (F^) is the gluon 
condensate and {qq) the light quark condensate. Taking the values as{F'^) ~ 
0.04 GeV^ pT], (gg) = (-270 MeV)^ [22)123] and A ~ 0.3 GeV its follows that 



4 



TZ ~ 6.8. Therefore, one expects 0" to be dominated by gluons, in agreement 
with our initial assumptions. 

The Effective Lagrangian . The effective theory includes the gluon field A^, 
the quark fields qf, / is a flavor index, an effective scalar field (6° belonging 
to the adjoint representation of SU (3) color group. The effective Lagrangian 
reads 

^ f 

+ ^ {D^r) {D^r) - Voctirr) + i^gf + + (^mw (4) 

where = dn + igT^A^ is the covariant derivative, stands for the genera- 
tors of SU (3) color group, m/ the current quark mass associated with flavor /, 
Voct is the effective potential associated with 0". Cqf is the gauge fixing part 
of the Lagrangian, Cgh contains the ghost terms and Cnew includes the new in- 
teractions. All terms in C are gauge invariant. The exception being Cqf- The 
coupling constant g is not the strong coupling constant. The effective model 
assumes that the non-perturbative physics is summarized in the color octet 
scalar field 0". Therefore, it is expected that the dominant contribution to 
any process to be associated with 0". The effective gauge coupling constant g 
parameterizes residual interactions and it should be a small number, allowing 
for a perturbative expansion in powers of g. 

The operators in Cmw can be classified accordingly to its mass dimensions. 
Given that 0" has dimensions of mass, the only gauge invariant dimension 
four operator that can be in Cnew without breaking parity is 



f 



r , (5) 



for a color octet scalar field 0". For a pseudoscalar 0", should be replaced by 

t"75. Note that we have considered a single coupling constant G4 to preserve 
flavor symmetry. The allowed dimension five operators in Cnew are 



G5E 

/ 



qq 



qf'q 



f 



q t"'j'^q 



f 



qf^rD^q 



(6) 



with 6*5 = g^/A and Fi = fi/A for i = 1,2,3 and where g^ and /j are di- 
mensionless coupling constants. For a pseudoscalar color octet field only the 
operators associated with coupling constant Fi and F2 should be modified by 
the introduction of a 75 between the quark fields. 
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The effective field theory is summarized in the Lagrangian (jl]), with C-new 
given by the sum of ([5]) and ([6]). £ includes the QCD Lagrangian and veri- 
fies the usual soft-pion theorems of chiral symmetry at low energy. The new 
interactions in Lnew introduce new vertices not present in the original QCD 
Lagrangian , namely qq(f)°'(f)^ and qq(f)°'A^^, which contribute to quark pro- 
cesses. Note that the only quark color singlet operator mimics the ^Pq model 
used to describe OZI-allowed mesonic strong decays. 

So far we have not yet specified the dynamics of the color octet field, i.e. have 
not defined its potential energy Voct- However, given that 0" cannot give rise 
to asymptotic S'-matrix states, for physical processes one has to integrate over 
0" and it is unnecessary to detailed its dynamics. As discussed bellow, it will 
be enough to specify its vacuum expectation values. 

Dynamical Gluon and Quark masses. The kinetic term associated with tp"" 
gives rise to a gluon mass term through the operator 

]^g^<P%T'^T\d<P''AlA'^^. (7) 

If ip"" acquire a vacuum expectation value without breaking color symmetry, 
i.e. 

(0'*) = and = , (8) 

and given that for the adjoint representation tr {t°'T^^ = Nc 5"*, then one can 
write the square of the gluon mass as 

ml = N,gW (9) 

where = 3. Note that the condensate {(f)°'(j)^), i.e. f^, and therefore the 
gluon mass is gauge invariant. The proof of gauge invariance follows directly 
from the transformations properties of 0". The gluon mass ([9]) is proportional 
to the effective gauge coupling g. As discussed previously, g should be a small 
number. However, this does not implies necessarily that is also small. The 
precise value of nig depends on the relative values of v and g. 



The same mechanism shifts the quark masses due to the operator 



qq 



(10) 



giving rise to a constituent quark mass 



l\j2 1 

= - (AT^^-l) G,v^ = m, - (H) 

For light quarks, the constituent mass is given by the quark self energy which, 
in the model, is linked with the gluon mass, the second term in equation (ITT]) : 
note, for our definitions, G5 is a negative number. If the constituent mass 
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o = + + ••• 

Fig. 1. Dressed quark propagator. The dashed hnes are associated with the color 
octet scalar field 4>'^. 

for the light quarks vanishes, chiral symmetry is recovered, and the model 
predicts that the gluon mass should also vanish. In this way, the model links 
chiral symmetry breaking with an effective gluon mass. 

Light Quark Condensate. The ligth quark condensate (qq) is an order param- 
eter for chiral symmetry breaking, whose computation requires the knowledge 
of the quark propagator. In lowest order, the quark propagator is given by the 
Feynman graph of figure [T], which includes the contribution from the conden- 
sate (0*^0''). The dressed quark propagator reads 

where M is the constituent quark mass, and the quark condensate is given by 



where Nf is the number of degenerate light flavors, uJ is an energy cut-off and 



J^{x,y) = xHn^^ +y\ (14) 

_|_ yZ 

For small M or for M <^uj, equation (|T3l) becomes 

(gg) = - ^^Muj\ (15) 

The quark condensate is proportional to the constituent quark mass, i.e. to the 
gluon mass squared, and to the square of the energy cut-off. A non- vanishing 
gluon mass means also a non-vanishing quark condensate and, therefore, chiral 
symmetry is broken. 

Equations (1151) and M when combined with the results of lattice QCD simula- 
tions [lOj allow us to estimate of the various parameters of the theory. Indeed, 
taking (qq) = (-270 MeV)^ M = 330 MeV it follows from that U = 626 
MeV for Nf = 2. For example, this is a typical value used for the cut-off 
when working with Nambu-Jona-Lasinio type of models. Taking the value for 
the gluon mass from Landau gauge lattice simulations rrig = 648 MeV, then 
G^/g^ = -0.295 GeV"^ and equation gives gv = 37A MeV. Note that gv 
is, apart the color factor y/Nc, the effective gluon mass and is of order Aqcd- 
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Testing the Gluon and Quark Mass relation. The effective model relates the 
constituent quark mass M and the gluon mass rUg through equation (fTT]) . 
For a vanishing current mass, equation f ill I) predicts a constant value for the 
ratio M/m^. This prediction can be tested looking at the solutions of the 
Schwinger-Dyson equations. In the following we will use the results published 
in |8]. For the gluon and ghost propagators, the authors use the results of 
lattice QCD simulations and solved the gap equation for a massless fermion. 
The calculation does not take into account fermion loops and can be viewed 
as a quenched approximation. 



For the gluon propagator, the inverse lattice propagator was fitted to 

' + pim?{q^)^ 



1 + 



96 7r2 



where 



2/ 2 

m [q 



m 



(16) 



(17) 



and m = 520 MeV, gf = 5.68, pi = 8.55 and p2 = 1.91 are the fitted parame- 
ters, /i = 4.3 GeV is the renormalization point and Ca = 3 is a SU{3) Casimir 
invariant. The gluon mass is given by the pole in the propagator and we will 
use the following definition 



rn^^q^ 



1 + 



96 



f,2 



(18) 



The fermionic gap equation was solved for two different ansatz for the quark- 
gluon vertex, namely a non-Abelian improved version of the Ball-Chiu vertex 
and an improved version of the Curtis-Pennigton vertex. The choice of vertex 
leads to slightly different quark mass. In order to distinguish, the results of 
the Ball-Chiu vertex will be referred as " EC" , while the results from using the 
Curtis-Pennington vertex will be referred as " CP" . 

The quark masses computed using the different vertex ansatz and the squared 
gluon mass are reported in figure [2l For a cut-off of = 0.626 GeV, i.e. up to 
momenta = 0.39 GeV^, the M/rrig ratio increases slightly. Indeed, M/m'^ 
changes by less than 15%, relative to its zero momentum value, for the BC 
vertex, and just below 19% for the CP vertex. 

We call the reader attention that the gluon and ghost propagators are quenched 
propagators and the quark mass was computed solving the fermionic gap equa- 
tion using the quenched propagators. A full dynamical calculation can change 
the above picture. Lattice QCD simulations [MpS] show a decreasing con- 
stituent quark mass M(g^) with the current mass in the infrared region. On 
the other hand, for dynamical simulations, the gluon propagator is also sup- 
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Fig. 2. On the left hand side, the plot shows the quark masses from solving the 
fermionic SDE gap equation, using different ansatz for the quark-gluon vertex, and 
the squared gluon mass computed from quenched lattice simulations. Note that 
M(g^) depends slightly on the definition of the quark-gluon vertex. On the right 
hand side, the plot shows the ratio M/rrig. 

pressed relative to the quenched propagator [25)126] . This suggests that rrig is 
enhanced when we approach the chiral limit. However, the dynamical ghost 
propagator is enhanced relative to the quenched propagator [25] in the infrared 
region. An enhancement in the ghost propagator implies an enhancement of 
the quark-gluon vertex [8j whose effects remain to be investigated. Despite 
the uncertainties related with the simulation of full QCD, hopefully, the sce- 
nario observed in figure Ej where the mass ratio changes by less than 19%, 
will remain essentially unchanged. Therefore, one can claim that the results 
summarized in figure [2] suggests that M/m^ is essentially constant in the non- 
perturbative regime. Certainly, given the uncertainties on the calculation of 
M and nig, figure [2] will not exclude such a possibility. 

Summary. We have explored the construction of an effective field theory de- 
scribing the low energy properties of QCD using scalar like fields as composite 
gluon operators. We found that color singlet scalar-like fields are unlikely to 
play a dominant role for the infrared properties. On the other hand, color 
octet scalar fields allow the building of SU{3) color gauge invariant effective 
interaction with interesting properties. The effective model is gauge invariant, 
flavor symmetric and satisfies the usual current algebra relations given that it 
includes the QCD Lagrangian . It assumes that the non-perturbative physics 
is resumed in the new condensate = v^6'^^, with f 7^ in the confined 

phase and t> = in the deconfined phase. 

Gauge invariance and the scalar condensate are enough to generate a gauge 
invariant gluon mass, which we associated with the transverse mass measured 
in Landau gauge Lattice QCD simulations. shifts the quark masses, jus- 

tifying the difference between current and constituent masses, and generates 
a light quark condensate proportional to the constituent quark mass. Further, 
the model predicts a constant M/m^ ratio. A massive infrared gluon prop- 
agator implies a non-vanishing constituent quark mass and a non-vanishing 
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light quark condensate, connecting with the breaking of chiral symmetry. 
Further, given the relation between nig and chiral symmetry, the model pre- 
dicts that the deconfinement transition and the restoration of chiral symmetry 
occur simultaneously. This is in good agreement with the predictions from the 
studies of the QCD phase-diagram. If the mass relations are unchanged by 
temperature and baryonic matter, the model provides a relation between the 
critical exponents for the gluon mass rjg and the light quark condensate rig 
near the chiral phase transition, namely rig = 2rig, that can be tested, for ex- 
ample, with lattice QCD simulations. The mass relation between M and rUg 
was tested against the decoupling solution of the SDE. 

The model includes a number of new interactions. The relations discussed 
above are the result of exploring two operators, the kinetic term for the scalar 
field and the only new color singlet operator at the quark level qq This 
last operator mimics the ^Pq model used to describe OZI-allowed mesonic 
strong decays. All the other operators seem to give marginal contributions to 
the strong decays because either (0") = 0, for the operators multiplying G4, 
F2 and F3, or dabb = 0, for the operator multiplying Fi. The new vertices can 
contribute the S'-matrix at the quark level and can give rise to new quark 
interactions. 

In QCD, the quark-gluon vertex is orthogonal to the gluon momenta. The 
vertices associated with the new interactions do not verify this property and 
can be viewed as mimicking the longitudinal part of the quark-gluon vertex, 
which is not constraint by Slavnov- Taylor identities. The effective model has 
phenomenological implications which have not been explored here and will be 
the subject of further publications. For example, the dressing of the quark- 
photon vertex due to (j)"' can give rise to quark anomalous magnetic moments 
of the type discussed in [28] . 

The model building blocks are the color octet field (f)"', the quark and the 
gluon fields. The hadronic matter is essentially quark matter. The experimen- 
tal situation for pure glue states is not so clear. One has to understand why 
there are no physical states than can be associated with 0". An explanation 
could be its low mass and the weak couplings to the quarks. Looking at the 
Lagrangian, one can estimate the contribution of the new terms to the mass of 
(j)" as 5ml ~ {qq)\- For G5 = -0.295 GeV'^ and {qq) = (-270 MeV)^, 
it follows that 5m^ < 76 MeV. Curiously, recently the authors of [29] have 
analyzed the BABAR bottomonium decay data and found an excess of signal 
which could be explained by a light scalar particle with a mass of 38 MeV. 

Finally, we note that the model, being a valid description for low energy QCD, 
predicts a behavior for the infrared gluon propagator which distinguish be- 
tween the so-called scaling solutions and decoupling solutions of SDE. The 
model is in close agreement with the decoupling like solution. 
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